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Abstract. 

We analyze nonequilibrium lattice models with up-down symmetry and two 
absorbing states by mean-field approximations and numerical simulations in two 
and three dimensions. The phase diagram displays three phases: paramagnetic, 
ferromagnetic and absorbing. The transition line between the first two phases belongs 
to the Ising universality class and between the last two, to the direct percolation 
universality class. The two lines meet at the point describing the voter model and the 
size ^ of the ferromagnetic phase vanishes with the distance e to the voter point as 
£ ^ e, with possible logarithm corrections in two dimensions. 
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1. Introduction 

Studies of phase transitions and critical phenomena in equilibrium and nonequilibrium 
systems have shown that the critical behavior can be organized into universality classes, 
which are identified by a small number of features among them the symmetry. The 
most widely studied universality class, due to its experimental implications, is the Ising 
universality class which includes equilibrium as well as nonequilibrium systems BE]- 
The main feature of the Ising universality class is the up-down symmetry. Another class, 
that includes only nonequilibrium systems, is the direct percolation (DP) universality 
class PHE]. The main feature that distinguishes this class from the others is the 
presence of a single absorbing state. That is, if a system displays a continuous phase 
transition to a single absorbing state with no extra symmetry or conserving laws then 
its critical behavior belongs to the DP universality class 

Our purpose here is to analyze models that possess features pertaining to the 
two universality classes. Here, we focus on models belonging to the Ising and DP 
universality classes, that is, systems with up-down symmetry and an absorbing state 
p El [ini dn [i2]. Due to the up-down symmetry, these systems, in fact, have two 
symmetric single absorbing states. Lattice models of this type have been studied 
numerically P revealing the existence of two transition points when a parameter is 
varied: a symmetry breaking point belonging to the Ising universality class and another 
belonging to the DP universality class. Another numerical study in two dimensions by 
the method of Langevin equation P, with two parameters, has shown that the two 
critical points evolve into two transition lines that meet at a certain point, given rise to 
a third transition line. This study conhrm the Ising and DP character of the two lines 
and the existence of a point on the third line corresponding to the voter model P, which 
is not the meeting point. It has been claimed that the symmetry breaking transition is 
not of the Ising type but Monte Carlo simulations on a two-dimensional interacting 
monomers model [12] conhrmed the Ising transition. In the present study, we wish to 
clear up this point and analyze the role played by the voter model. This is achieved 
by studying lattice models with up-down symmetry and two absorbing states in two 
and three dimensions by numerical simulations and mean field approach. Our results 
conhrmed the existence of the Ising and DP lines and revealed that two lines meet at 
the voter point and that this is the only point at the third line which is critical. 

The phase diagram in the two-dimensional space of parameters displays three 
phases: paramagnetic (P), ferromagnetic (F) and absorbing (A), as shown in £gure[T^. 
The PF line that separates the P and F phases is a critical line belonging to the Ising 
universality class. The FA line that separates the A phase and the F phase, understood 
as an active phase, is a critical line belonging to the DP universality class. The two 
lines, PF and FA, meet at the point corresponding to the voter model. As to the PA 
line, that separates the P and A phases, it is a discontinuous transition line that ends 
at the voter point. Along the PA line the jump in the order parameter is a constant 
and equals the jump observed in the voter model. For an appropriate description of the 
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(a) 




Figure 1. (a) Phase diagram in the field parameters e and r. The full circle, located 
at the origin, corresponds to the voter model (V). The phases are: paramagnetic (P), 
ferromagnetic (F) and absorbing (A). The solid lines are continuous phase transitions 
belonging to the Ising and DP universality classes and the dashed line is a discontinuous 
phase transition, (b) Magnetization m as a function of e and r. 


critical properties around the voter point, it is convenient to introduce two parameters, 
£ and r, that dehne a reference frame centered with origin at the voter point. The £-axis 
and r-axis are parallel and perpendicular to the direction dehned by the discontinuous 
transition line, as shown in hgure [1^. 

The order parameter m of the F phase, shown in hgure [Tb is assumed to obey the 
following scaling relation, aronnd the voter point, 

m =(1) 

where 0 is the crossover exponent and '0+ and ip- are valid for e > 0 and £ < 0, 
respectively. The nniversal fnnction 'ip-{x) is a step fnnction. The nniversal fnnction 
'ip+i^x) vanishes according to t/’+(3^) ~ (a — where (di is the Ising order parameter 
exponent and a is a positive constant so that, near the Ising critical line, the order Ising 
parameter m behaves as 

m = A{ei-ef\ (2) 

where Si = ar^^d ^nd the amplitnde A diverges as A ~ universal fnnction 

'ip+{x) approaches its maximum value 1 according to 1 — V’+ ~ where /Sdp is 

the DP order parameter exponent and 6 is a constant so that, near the DP critical line, 
the DP order parameter p = 1 — m behaves as 

p = 5(5 + £Dpf°^ (3) 

where £dp = br^ and the amplitnde B diverges as i? ~ nf-hYyvl<t>_ 

The Ising and DP critical lines occur when m —)■ 0 and p —)■ 0, respectively, and 
are given by £ = ar^ld and e = —br^ld so that the size of the F phase behave near the 
voter point as 




( 4 ) 








Nonequilibrium lattice models with two symmetric absorbing states 


4 


Table 1. Probability of transition ai ——tTi. 
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2. Model 


We consider stochastic variables that take two values, cTj = ±1, called spin variables. 
They are located on the sites of a regular lattice which is either a two-dimensional 
triangular lattice or a three-dimensional cubic lattice. Notice that both lattices have 
the same coordination number, which is six. The system evolves in time according to 
a continuous time Markov process with a single-site change. The single-site transition 
rate Wi{a) from a = (ai,..., cTj,..., ctat) to a* = (di,..., —cxj,..., ctat) is set up in such 
way as to hold the properties: a) it has the up-down symmetry, that is, it is invariant 
under the transformation a —)■ —a; b) it depends on the neighboring spin variables 
only through the sum of these spin variables. A transition rate that fulhls these two 
properties is given by 

- 0-i/(Si)}, Si = Y, (^i+5, (5) 

where the sum in S extends over the six nearest neighbor sites, and /(s) is an odd 
function of s. 

The most general form of f{s) has three parameters. Let us denote by aj, 
j = 1,2, 3,4, 5, 6 the spins of the six nearest neighbor sites of a central site j = 0, 
so that 


So — C’’l+C’’2 + C’'3 + C’'4 + C’’5 + C’’6- 


( 6 ) 


Taking into account that aj = ±1, then sq may take only the values 0, ±2, ±4, ±6 from 
which follows that the most general form of an odd function /(s) can be written as 

f(s) = -s + -s^-\ -s^, (7) 

where a, b and c are parameters. Here, it is convenient to dehne new parameters 
Po = [1 + /(6)]/2, pi = [l + /(4)]/2 and p 2 = [1 + /(2)]/2, that is. 


Po — 2 243c), 

Pi = —(1 -|- 2a -|- 86 -|- 32c), 


( 8 ) 

(9) 


P 2 = -{l + a + b + c). 

In terms of these parameters, the transition probabilities are as shown in table [H 


( 10 ) 
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Two well known models are special cases. When po = 1) Pi = 5/6 and p 2 = 2/3, 
(a = 2/3, b = c = 0) we recover the so called voter model [21 [9] whose transition rate is 
given by 

u;™ter(^) ^ g. = (11) 

A 6 g 

The Glauber-Ising model, whose transition rate is given by 

tanh(i^Si)}, Si = Y^ ai+s, (12) 

^ <5 

is recovered when the parameter are connected by 

Po =-(1 + lanhGit'), (13) 

Pi = ^(1 Ttanhdit'), (14) 

P 2 = ^(1 + tanh2it'). (15) 


In this case, detailed balance is obeyed and the stationary distribution is the Gibbs 
distribution associated to the Ising Hamiltonian. In other cases there is no detailed 
balance as can be seen by the following sequence of transitions. Gonsider two nearest 
neighbor sites initially in states (++). All the other nearest neighbors of this pair of 
spins are positive except one of them which is negative. Next, they change their states 

according to the cyclic sequence (++) (H—) (-) (—h) —>■ (++)• Using the 

rules of table [T] we find that the ratio between the transition probabilities of this path 
and of its reverse equals (1 — Pi)^P 2 Po/(l “ Po)(l ~ P 2 )PiPi, which in general is distinct 
from the unit. Thus, there is no microscopic reversibility and no detailed balance. 

To meet the two features of the Ising and DP universality class, that is, up-down 
symmetry and two equivalent single absorbing states, the transition rate should lead to 
an absorbing state. This is accomplished by imposing the restriction pq = 1 because in 
so doing the transition rate will vanish whenever all spins are up or all spins are down, 
which are the two states identihed as the equivalent absorbing states. In this case, 
the transition rates do not obey detailed balance, which amounts to say that the two 
features mentioned above are only possible in systems that in the stationary states are 
out of equilibrium. When po = 1) fhe phase diagram is restricted to the plane (pi,P 2 ), 
and the voter model corresponds to the point (5/6, 2/3) of this diagram. As we shall 
see, the phase diagram displays a paramagnetic phase, a ferromagnetic phase and an 
absorbing state, as explained before. 

3. Mean-field approximation 

A qualitative phase diagram is obtained by the use of a mean-held approximation. To 
this end, we start by writing the function /(s) in the form 

, , A B C 

J\^) T T 


( 16 ) 
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Figure 2. Phase diagram in the plane {pi,P 2 ) according to mean-field approximation. 
The full circle, located at {pi,P 2 ) = (5/6, 2/3), corresponds to the voter model. The 
phases are: paramagnetic (P), ferromagnetic (F) and absorbing (A). The solid lines are 
continuous phase transitions and the dashed line is a discontinuous phase transition. 


where the hrst summation runs over i from 1 to 6, and has 6 terms; the second 
summation runs over i,j, k with the restriction l<i<j<k<6 and has 20 terms; and 
the third summation runs over i, j, k, i, n with the restriction l<i<j<k<^<n<Q 
and has 6 terms. The parameters A, 5, and C are related to po; Vi) and p 2 by 

Pq =-{1 + A + B + C) (17) 

P. = 1(1 + - \c) (18) 

(19) 


The time evolution of the magnetization (uo) is given by 

^Mo) = -2(uoWo(a)) = -(uo) + (/(so)) (20) 

dt 

or 

^(uo) = -(uo) + A{ai) + B{aiajak) + Ciaiajakaian) (21) 

where translation invariance has been used. Next, we use the simplest mean-field 
approximation, for which {a^ajak) = and {aiajakCJ^an) = m®, where (uj) = m, 
to write 

= —(1 — A)m + Bm^ + Cm^. (22) 

CL L 

Here, we are interested in the case po = 1 which is equivalent tol — H = i?-|-C'so 

that 

uTD 

-—=—{B + C)m + Bm^ + Cm^, (23) 

CL L 
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which is identical to the deterministic part of the Langevin equation considered in 
reference [9]. In the stationary state, the possible solutions are as follows, a) m = 0, 
corresponding to the paramagnetic (P) phase; it is stable as long as S + C > 0, which is 
equivalent to 12pi + 15p2 < 20. b) m = ± 1 , corresponding the absorbing (A) phase; it 
is stable as long as B + 2C < 0, which is equivalent to pi > 5/6. c) m 7 ^ 0 and m 7 ^ ±1, 
corresponding to the ferromagnetic (F) phase. The phase diagram in the plane {pi,P 2 ) 
is shown in figure |2l The P-F line is given by i? + C = 0, or 12pi + 15p2 = 20, and 
the F-A line is given by i? + 2C = 0, or pi = 5/6. The point corresponding to the 
voter model is located at the point where the two lines meet, that is, at i? = C = 0 , or 
Pi = 5/6 and p 2 = 2/3. The P-F and F-A transitions are continuous whereas the P-A 
transition is discontinuous. 

The discontinuous P-A transition line shown in figure [2] was obtained as follows. We 
start by writing equation fl2^ in the form dm/dt = —df /dm where / = {B + C) 12 — 
Bm^/A — Cm^/Q -|- K. Interpreting /(m) as a free energy then the discontinuous 
transition occurs when /(O) = /(±1), which gives ?>B -|- AC = 0 or 5p2 + 8pi = 10. 

Inside the F phase, the magnetization m is given by 

1/2 


m = 


( Ibp2 + I2pi-2Q \ 






lbp2 - I2pi 

Near the critical line PF, pi = ppp = (20 — 15p2)/12, we may write 


m = A{pi 


A = 


P2 - 


1 - 1/2 


(24) 


(25) 


so that the amplitude A diverges with an exponent 1/2 as one approaches the voter 
point. Analogously, near the FA transition, pi = ppA = 5/6, the order parameter p 
behaves as 


P = B(Pfa -Pi), 



(26) 


SO that the amplitude B diverges with an exponent 1 as one approaches the voter point. 
It is worth to notice that the size £ = p-pA — Ppa of the ferromagnetic phase, given by 


£ = 



(27) 


vanishes linearly as one approaches the voter point. 


4. Simulations 

The results obtained from numerical simulations on the cubic and triangular lattice allow 
us to conclude that the topology of the phase diagram is the same as the one obtained 
from the mean-field approximation. The numerical simulations on these lattices were 
carried out as follows. At each time step a site of the lattice is chosen at random 
and updated according to the transition probabilities in table [H restricted to the case 
Po = 1. The averages of several quantities were obtained in the stationary state. Defining 
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Figure 3. (a) Magnetization m versus pi along P 2 = 1- The transition occurs at 
Pi = 0.5930 as can be inferred from the plot (b) of im}/^ versus pi, where j3 = 0.326. 


s = (1/iV) where N is the number of sites, then the quantities of interest to analyze 
the ferromagnetic-paramagnetic line are the magnetization 


m = (|s|), 
the susceptibility, 


(28) 

(29) 


the Binder cumulant 


U = 1- 


3(s2)2- 


(30) 


Examples of these quantities, for the cubic lattice with N = sites, calculated 
along the line p 2 = 1, are shown in hgures [3] and 01 At the P-F transition point, the 
Binder cumulant f/ is a constant, a result that allows to located the P-F transition line. 
The critical behavior of m and y allows to say that, with respect to the ferromagnetic- 
paramagnetic transition, the model belongs to the universality class of the Ising model. 

To analyze the ferromagnetic-absorbing transition, where the ferromagnetic phase 
is identihed as the active phase, we dehne the quantity p = 1 — m, which is the order 
parameter that characterizes the transition from active to absorbing. In hgure [U we 
show p versus pi along the line p 2 = 1- The location of the transition is obtained by the 
extrapolation p —)■ 0, what allows us to draw the F-A line. The critical behavior shows 
that with respect to this transition, the model belongs to direct percolation universality 
class. 

Using the data of the P-F and F-A transitions obtained for several values of p 2 we 
draw the phase diagram of hgure [6^. As expected the two lines, P-F and F-A, meet 
at the voter point located at pi = 5/6 and p 2 = 2/3. In hgure [6b we show the size 
Api = Pdp — Pi of the F phase as a function of the distance from the voter point 
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Figure 4. (a) Susceptibility x and (b) Binder cumulant U versus pi along p 2 = 1. 
The transition occurs at pi = 0.5933. 


Ap 2 = P 2 — 2/3. As can be seen in this figure, the data points has a finite slope at the 
origin which amounts do say that the scaling relation 

Api ~ (Apa)'^, (31) 

is fulfiled with 0 = 1.0. 

We have also determined the nature of the phase transition between the P and A 
phases. As we cross the transition line the order parameter jumps from zero, which is 
the value of m within the paramagnetic phase, to the maximum value m = 1, which 
is the value of m inside the absorbing phase. The susceptibility y increases as one 
approaches the transition line from the paramagnetic phase and reaches a finite value 
Xo{Ap 2 ) which depends on the distance Ap 2 = P 2 — 2/3 from the voter point. As one 
approaches the voter point xo increases without limits and diverges with an exponent 
7 = 1. Notice that, according to our numerical results for the susceptibility obtained 
along p 2 = 2/3, the susceptibility y diverges as one approaches the voter point from the 
paramagnetic phase with the same exponent 7 = 1 . 

We have also studied the model defined on a triangular lattice by using the same 
methods. As can be seen in figure [TJ the phase diagram also displays the three phases 
found in the cubic lattice. However, the ferromagnetic phase is very narrow, becoming 
very difficult to extract the crossover exponent but the two lines, PF and FA, seems to 
meet tangentially in consistency with an exponent 0=1 with logarithm corrections as 
we will argue below. 

5. Discussion 

We have studied the phase diagram of a system belonging to the Ising and DP 
universality classes. The phase diagram in the two parameters has an Ising critical line 
separating the paramagnetic and the ferromagnetic phase, a DP critical line separating 
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Figure 5. (a) Order parameter p of the F-A transition as a function of pi for p 2 = 1, 
The critical point occurs at pop = 0.7670 as can inferred from the plot (b) where 
h = 0.81 




Figure 6. (a) Phase diagram in the plane {pi,P 2 ) according to simulations on a cubic 
lattice. The phases are: paramagnetic (P), ferromagnetic (F) and absorbing (A). The 
solid lines are continuous phase transitions and the dashed line is a discontinuous phase 
transition. The three lines meet at the point {pi,P 2 ) = (5/6, 2/3), corresponding to 
the voter model, (b) Linear size Api = ppA — Ppf of the ferromagnetic phase as a 
function of the distance from the voter point Ap 2 = P 2 —2/3. The dashed straight line 
shows that the data points have a finite slope at the origin. 


the ferromagnetic phase and the absorbing phase, and a discontinuous phase transition 
line separating the paramagnetic and the absorbing phase. The three lines meet at the 
voter point in such as way that the size of the ferromagnetic phase vanishes as one 
approaches the voter point with a crossover exponent 0. The crossover exponent should 
be understood as the ratio between the two exponents related to the two relevant scaling 
helds e and r around the voter point. Above the critical dimension dc of the voter model, 
these two exponents do not depend on dimension so that the same should happen to 
the crossover exponent. Taking into account that the upper critical dimension of the 
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Figure 7. Phase diagram in the plane {pi,P 2 ) according to simulations on a triangular 
lattice. The full circle, located at (pi,_P 2 ) = (5/6, 2/3), corresponds to the voter model. 
The phases are: paramagnetic (P), ferromagnetic and absorbing (A). The inset shows 
an enlargement around showing the ferromagnetic phase between the P-F and F-A 
transition lines. 


voter critical point is dc = 2, we expect the crossover exponent 0 to have the same value 
for (i > 2 with possible logarithm corrections when d = 2. The results coming from 
mean-held theory gives 0=1 which allows to conclude that 0 = 1 for d > 2, a result 
conhrmed by our calculations in d = 3. As to the two-dimensional case, the size of the 
ferromagnetic phase is too narrow for a precise numerical calculation but is consistent 
with the result 0=1 with logarithm corrections. 

The results for the susceptibility around the Ising and DP lines are in accordance 
with these two universality classes. Concerning the discontinuous transition from the 
paramagnetic to the absorbing transition, we found that the susceptibility is hnite at 
the transition line but diverges as one approaches the voter point. 

It would be very interesting if the picture that we have drawn concerning the 
critical behavior of systems belonging to the Ising and DP universality classes could 
be conhrmed by a renormalization group calculation by means of a Langevin equation 
formulation. 
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